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Abstract 

Radioactive decay of an unstable isotope is widely believed to be exponential. This 
view is supported by experiments on rapidly decaying isotopes but is more difficult to 
verify for slowly decaying isotopes. The decay of 14 C can be calibrated over a period of 
12,550 years by comparing radiocarbon dates with dates obtained from dendrochronology. 
It is well known that this approach shows that radiocarbon dates of over 3,000 years are in 
error, which is generally attributed to past variation in atmospheric levels of 14 C. We note 
that predicted atmospheric variation (assuming exponential decay) does not agree with 
results from modelling, and that theoretical quantum mechanics does not predict exact 
exponential decay. We give mathematical arguments that non-exponential decay should 
be expected for slowly decaying isotopes and explore the consequences of non-exponential 
decay. We propose an experimental test of this prediction of non-exponential decay for 
14 C. If confirmed, a foundation stone of current dating methods will have been removed, 
requiring a radical reappraisal both of radioisotope dating methods and of currently pre- 
dicted dates obtained using these methods. 



1 Introduction 



Radioactive decay is almost universally believed to satisfy the exponential decay law over many 
half lives and in particular cases where it has been tested [TJ [2] it has been found to accurately 
represent the decay of an unstable isotope. However, as nuclear physics as a discipline is 
little over a century old, these experiments have necessarily been restricted to rapidly decaying 
isotopes. This principle of exponential decay is then assumed to hold in many other situations 
where there is no direct experimental verification, notably for slowly decaying isotopes which 
are widely used for dating objects ranging from a few hundred years old, to billions of years 
old. Newton summarised this situation as follows: 'Although in some instances the exponential 
law has been experimentally well verified over many lifetimes, it has certainly not been checked 
in many cases in which it is nevertheless assumed to be valid' [31 p608]. 

So is there any experimental or theoretical evidence that can be used to verify, or disprove, 
this widely accepted assumption for slowly decaying isotopes? 
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With regard to experimental evidence, we consider in detail the decay of 14 C which is the 
basis of radiocarbon dating. We propose that non-exponential decay may be a contributing 
factor in the known discrepancy in radiocarbon dates of over 3,000 years. We also consider 
the theoretical aspects of radioactive decay as described by quantum mechanics, and again 
present arguments for non-exponential decay for slowly decaying isotopes. We then consider 
the consequences if non-exponential decay is assumed and conclude with a proposal for an 
experimental test of non-exponential decay. 

2 Radiocarbon dating and calibration 

The radioisotope 14 C is the basis for radiocarbon dating in which it is assumed that living 
organisms have a 14 C/ 12 C ratio which is the same as that in the atmosphere and that this 
atmospheric ratio has been constant in the past. The interaction between the organism and the 
atmosphere stops at death, and the 14 C that remains in the organism decays. By measuring 
the 14 C/ 12 C ratio directly (using Accelerator Mass Spectrometry (AMS)) or the rate of decay 
of 14 C (using radiometric methods), the time since death can be estimated. Clearly there is 
another assumption underlying this process, namely that the decay of 14 C is exactly exponential 
for all time. The half-life for 14 C was initially found to be 5, 568 ± 30 yr (the Libby half-life) 
but was later changed to 5, 730 ± 40 yr (the Cambridge half life) [I]. The Libby half-life is still 
used by convention in calculating 14 C ages [U EJ [7J. 

Radiocarbon dates can be calibrated using dendrochronology in which the age of a sequence 
of trees of increasing age can be determined very accurately by counting the growth rings and 
cross-dating different samples of overlapping age [I]. There can be occasional problems with 
missing rings or double rings produced in a single year, but with careful analysis and checking 
of several samples the errors in this process are very small ('on the order of 1 yr' |8J). The 
longest continuous tree-ring chronology of 12, 593 years is based on German oak and pine trees, 
supplemented by pines from Switzerland [9j[T0]. This approach is unique in providing accurate, 
absolute dates over a period of many thousands of years. 

The dates of a tree sample obtained using dendrochronology and by radiocarbon dating 
are found to be very similar over the initial 3, 000 year period, but for older trees there is 
a significant difference between these two dates of up to 15%. It is assumed that the date 
obtained from dendrochronology is the more accurate one and so these two dating methods 
can be used to generate a calibration curve for converting radiocarbon dates to calendar dates. 
The currently accepted calibration curve is IntCal09 [11] which is based on dendrochronology 
for the first 12,550 years and extends back to 50,000 years using other methods. We note that 
radiocarbon dates have units BP (before present, where the present is the base year of 1950) or 
kBP, while calibrated dates obtained from the calibration curve have units cal BP or cal kBP. 

While the use of the calibration curve may have resolved the issue for the practitioner, who 
now has a valuable and accurate dating tool, for the theorist it leaves the obvious question as 
to what is the cause of this discrepancy in the radiocarbon dates? 

3 Atmospherics 

The discrepancy in the calibration curve is generally attributed to past atmospheric variation in 
the 14 C/ 12 C ratio, which contradicts the assumption made in the dating process that the ratio 
was constant in the past. Four factors are considered to affect this atmospheric ratio, namely 
primary cosmic ray flux, strength of the solar electromagnetic field, terrestrial magnetic field 
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intensity and the structure of the carbon cycle [U Q2J [33] . The combined effect of these factors 
over the last 12, 550 years is not known and since we do not have an independent method of 
measuring the atmospheric ratio over many thousands of years, it is not possible to correlate 
the variation in the atmospheric ratio with the discrepancy in the calibration curve. However, 
this is a significant assumption made in the carbon dating process which is known not to hold 
and so is generally considered to be the sole cause of the discrepancy in the radiocarbon dates. 

If the decay of 14 C is assumed to be exponential, then the discrepancy in the calibration 
curve can be used to predict the variation A 14 C of the atmospheric ratio in the past [TJ]. Data 
from a stalagmite [13] and from marine sediments [15] have been used to construct a record of 
A 14 C going back to 45 or 50 cal kBP respectively. These records show much larger variations 
than are found over the 12.55 cal kBP period that we have considered. Both studies also used 
various box models of the global carbon cycle, which typically did not achieve levels of A 14 C as 
high as those derived from the data unless extreme values of the parameters were used. Chiu 
et al. [16] analysed these models and stated that 'there is no commonly accepted explanation 
for the high atmospheric A 14 C values recorded in most archives'. They also concluded that 
'the discrepancies between measured A 14 C and modelled A 14 C remain unresolved'. They noted 
that one way to reduce the predicted high levels of A 14 C would be to increase the half life of 
14 C. A half life of 6,030 years gives results that are 'entirely consistent with the Beck et al. 
(2001) model (A) prediction' [16], which corresponds to an increase in the half life of over 5%. 
However, recent experimental results have proposed that the half life should be decreased by 
2±1% [T7], which would of course increase the predicted A 14 C values even further, and so there 
is no justification for the increase in the half life that is required to reduce the high A 14 C values. 

While the date range used by these authors is way beyond the range that we are considering, 
these studies suggest that variation in the atmospheric ratio may be insufficient on its own to 
explain the discrepancy in the calibration curve. However, allowing for non-exponential decay 
as well as atmospheric variation would mean that the predicted A 14 C curve could be reduced 
to be more in line with the predictions obtained from the modelling. So is it possible that some 
of the discrepancy in the radiocarbon dates could be due to non-exponential decay of 14 C? 

4 A statistical perspective 

Before considering the possibility of non-exponential decay, we first review the common simple 
derivation of exponential decay. 

Radioactive decay occurs when an unstable atomic nucleus spontaneously emits energy and 
matter, often transforming into a new element in the process. The unstable 14 C nucleus can 
undergo radioactive beta decay in which an electron and an anti-neutrino are emitted and the 
14 C decays to the stable isotope 14 N. This decay process is assumed to be stochastic at the 
atomic level and hence unpredictable, but the average rate of decay for a large number of 
atoms is often highly predictable. Assuming that each atom is equally likely to decay at any 
time and that the probability of decay is independent of the age of the atom leads to a Poisson 
process [18] which can be described in terms of the differential equation 

A = -aA, A(0) = A , (1) 

where A(t) is the 14 C/ 12 C ratio at time t with initial value Aq at time t — 0, which corresponds 
to the atmospheric ratio, a is the decay constant which determines the rate of decay and 
A = dA/dt. The solution of equation is 

A(t) = A e- at , 
3 



which gives the expected exponential decay. 

We note that this derivation of exponential decay is based on a simple statistical argument 
and makes no attempt to understand the mechanism that causes the decay. It has been noted 
that 'the exponential law. . . is experimentally verified in some cases to an astonishing degree, 
much more so than one has a right to expect considering the rather flimsy basis on which it is 
usually derived' [31 p594] . So what is the correct theoretical basis for understanding radioactive 
decay? 

5 A quantum mechanical perspective 

The theory of quantum mechanics describes the process of radioactive decay in terms of solutions 
of the Schrodinger equation [T9J [20] . This theory highlights one of the deficiencies of the simple 
equation (pQ), namely that it does not contain history-dependent terms, as are contained in the 
equation for the probability amplitude derived from the theory of quantum mechanics [191 E] • 

The general solution of the Schrodinger equation using first order perturbation theory does 
not give the probability of non-decay as a simple exponential. However, if some approximations 
are used, known as the Fermi Golden Rule [191 120] , then this probability amplitude becomes 
exponential. Thus, 'the exponential decay law, for which we have so much empirical support 
in radioactive decay processes, is not a rigorous consequence of quantum mechanics, but the 
result of somewhat delicate approximations' [201 p513]. 

The theory of quantum mechanics predicts variations from pure exponential decay at very 
short times and at very long times [21] . However, these deviations occur in regimes that are 
currently beyond the reach of experimental verification [2j. 

Thus, there is no theoretical basis from quantum mechanics for exponential decay in all 
circumstances. We might therefore ask why the theoretical possibility of non-exponential decay 
has not been observed experimentally. 

6 A mathematical perspective 

One important mathematical property of the differential equation ([T]) is that it is linear. Ian 
Stewart observed that 'today's science shows that nature is relentlessly nonlinear' [251 p83]. 
This does not mean that there is no room for linear equations but if the world is essentially 
nonlinear, then linearity is a special case and as such must be justified for a particular problem. 
So could this linear equation in reality be an approximation of a nonlinear equation? And what 
difference would it make if this equation was actually nonlinear? 

To address these questions, we note that in the context of radioactive decay, there are two 
quantities that can be used to distinguish four scenarios namely (i) the rate of decay and (ii) 
the timescale (relative to the half life) over which measurements are made. This is illustrated 
in fig. [TJ We discuss the terms "rapid decay" and "slow decay" in more detail later. 

We first consider the short time case. Plotting A(t) on a log scale results in a straight line 
graph if A(t) is exponential. On a timescale that is small relative to the half life, only a short 
section of the line is observed and this may be approximately straight, even if it is not straight 
over a much longer timescale, so that decay, when considered over a short time period, appears 
to be exponential. This covers the two possibilities in the bottom row of fig. [H 

In quantum mechanics, the decay of an unstable state is described by the Schrodinger 
equation, from which a convolution differential equation can be derived for a complex transition 
amplitude z(t) from the initial unstable state to itself. The probability of finding the system 
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Figure 1: The type of decay that is likely in four different cases. 

still in the initial state at time t is |z(t)| 2 , which is the quantity that is usually assumed to 
decay exponentially [20j 12] • 

The Schrodinger equation is linear, which follows from the Principle of Superposition [20J, 
as is the derived convolution differential equation for z(t). However, if we work with the two 
real variables A = \z\ 2 and 9 = axg(z), then the differential equations in these two variables 
are nonlinear. Scaling A by 1/Aq ensures that it satisfies the initial condition A(0) = 1. The 
equation for A can then be expressed in the form 

A = -aA + g(A,9,t), (2) 

where the nonlinear function g involves an integral over past time and satisfies g(0,9,t) = 0. 
We are not concerned with the equation for the second variable 9. 

We now define rapid decay to be when a is large compared to the nonlinear term g(A, 9, t) 
and slow decay to be when the two terms are of comparable magnitude. 

1. If decay is rapid, then the nonlinear term in equation (J2J) will be small compared to the 
linear term resulting in decay that is very close to exponential for a long time. This covers 
the top left box of fig. [TJ 

2. If decay is slow, then the linear and nonlinear terms will be more comparable in magnitude, 
and the nonlinear term will then play an important part in the decay dynamics, resulting 
in non-exponential decay. This is the remaining combination in fig. [TJ 

The above discussion contains a number of assumptions which cannot easily be verified. 
However, it is undoubtedly true that accurate measurements have been made in three cases, 
as also indicated in fig. [TJ with all such measurements strongly supporting the assumption of 
exponential decay. For example, the decay of 56 Mn, which has a half life of 2.5785 h, has 
been measured experimentally up to 45 half lives in order to search for non-exponential effects, 
but none were found (2]. In the fourth case, it has not been possible to accurately measure 
slow decay over long time periods, which is precisely the case where it is hard to justify the 
assumption of a linear equation giving rise to exact exponential decay. 

Thus, it is possible that the slow decay of 14 C could be non-exponential, and consequently 
that this non-exponential decay is a contributing factor to the discrepancy in the radiocarbon 
dates described above. 
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Figure 2: The calibration curves Cp (red) and Cams (blue) together with the line r r = {acl ( ^L) T d 
(magenta). 

7 What if? 

We have thus far made a case that variation in the atmospheric ratio may be insufficient on 
its own to explain the discrepancy in the calibration curve, and that the possibility of non- 
exponential decay should also be considered. We therefore now assume that the discrepancy 
in the calibration curve is a combination of atmospheric variation and non- exponential decay. 
Of course our big problem, having made this assumption, is that we do not know how much 
each of these two factors contributes to the discrepancy. There are two extremes that could be 
considered, namely (i) that the discrepancy is entirely due to atmospheric variation (and hence 
that decay is exponential) or (ii) that there is no atmospheric variation and the discrepancy 
is entirely due to non-exponential decay. Case (i) is usually assumed to apply, and the conse- 
quences of this assumption have been studied, resulting in the A 14 C curve which we discussed 
above. We now make the alternative extreme assumption described by case (ii), in order to 
consider the consequences of this assumption, whilst recognising that this does not occur in 
practice. 

Before continuing, we digress briefly to consider the INTCAL09 dataset [11] which consists 
of an underlying curve with superimposed small oscillations, known as "wiggles". A close 
correspondence between these wiggles and solar activity related to sunspots over a period of 
400 years has been observed [T2J [22] ■ Thus, we ignore the wiggles and extract the underlying 
trend in the data. Let r r and r d be the dates of a tree sample obtained using radiocarbon dating 
and dendrochronology respectively. We fit the data with a polynomial of degree 12 which we 
denote by r r = Cpijd). This smoothed curve is shown in fig. [2j 

There are two methods used for carbon dating, namely radiometric methods (/3-counting), 
which measure the rate of decay, and AMS, which measures the 14 C/ 12 C ratio directly. A 
consequence of the assumption in case (ii) is that these two methods will give different predicted 
dates for a given sample. To quantify this difference, let Y(t) be the (non-exponential) decay 
profile of the 14 C/ 12 C ratio which satisfies Y(0) = A , where A is the (assumed constant) value 
of the atmospheric ratio. We also define a non-dimensional decay function y(t) = Y(t)/Ao which 
is the proportion of the original 14 C that remains at time t, and satisfies the initial condition 




1/(0) = 1. 

With radiometric methods, the sample activity, or rate of decay, is measured, normalised to 
take account of isotopic fractionation, and adjusted to the 1950 level. This normalised sample 
activity, denoted by A^, is compared with the similarly normalised activity of an oxalic acid 
standard, denoted by A ON , which is considered to be the same as the activity of the atmospheric 
carbon. The radiocarbon age r r (BP) of the sample is then given by [5] 

r r = -i-ln(4si), (3) 

OIL VAon/ 

where oll — In 2/5568 yr -1 , which involves the Libby half life. 

If the decay curve is given by Y(t), then the activity of the sample at time to is v4sn = —Y{to)- 
The activity of the standard, assuming that decay is exponential initially, is Aqn = olcAq, 
where ac = ln2/5730 yr^ 1 , which is the decay constant based on the more accurate Cambridge 
half-life. Substituting for A SN and A ON in ([3]) then gives 





(4) 



Rearranging equation (j4j) gives 



y(t ) = -a c e~ 



■a L T r 



The calibration curve IntCal09 in the range of interest of 0-12.55 cal kBP was derived almost 
exclusively using radiometric methods. Using the smoothed calibration curve we therefore have 
that r r = Cp{j d ) where r d is the true age of the sample obtained by dendrochronology. This is 
the same as the time to, and so we have 



y(t ) = -a c e 



-a L Cp{t ) 



This equation holds for all values of to, and so we replace it with t. Integrating and using the 
initial condition y(0) = 1 then gives 

y (t) = 1 - ac t e' a ^ (s) ds. (5) 

JO 

This function is compared with the corresponding exponential decay curve e~ act in fig. [31 from 
which it can be observed that there is very good agreement between the two curves initially 
with more of a difference developing after approximately 6,000 years. 

The ratio of the function y(t) and the exponential e~ act shows small amplitude oscillations 
initially (see fig. H|). These arise due to the fact that the calibration curve initially oscillates 
also (see fig. [2]). In a simple tunnelling problem, the survival probability was found to be non- 
exponential, and also showed oscillations [6], although the amplitude was much larger than 
those shown in fig. |H 

We note that if it is assumed that decay is exponential, then substituting y(t ) = e~ act ° 
into (jlj) and replacing to with the dendrochronologically determined age T d gives 

r r = —r d = 0.9717r d , 

which is slightly lower than the line r r = r d . The reason for this is the historic use of rather 
than ac in (El). 
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a c = In 2/5730 yr _1 . 
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Figure 4: Plot of y(t)/e act showing small amplitude oscillations. 

Using AMS, both the 14 C/ 13 C and 13 C/ 12 C ratios of the sample and a standard are measured. 
The 14 C/ 13 C ratio is normalised to take account of isotopic fractionation, which adjusts the 
13 C/ 12 C ratio of the sample to be the same as that of the standard, and an adjustment to the 
1950 level is made. The normalised 14 C/ 13 C ratios of the sample and the standard are denoted 
by (14/13)sn and (14/13)on respectively and the radiocarbon age r r (BP) is then given by [7] 

1 /(14/13)sn\ 
a L \(14/13)on/ 

where, by convention and for consistency with (j3J), the Libby half life is again used for the 
decay constant. Again assuming that the true decay curve for the 14 C/ 12 C ratio is given by 
Y(t), then we have that (14/13)sn/(14/13)on = Y(t )/A , where t is the true age, since the 



s 



C/ 12 C ratios of the sample and the standard are the same. Thus, equation (E]) becomes 



T r = In 



( 



Y{to) 



A 



) 



lay (t ). 



(7) 



Now to is the true age of the sample, which we denote by to — Tj, and so from equation ([7]) we 
now get a new calibration curve based on AMS measurements given by 



where the function y is given by (J5J). This curve is shown in fig. [H together with the line 
Tr = (oic/otLjTd which is associated with exponential decay. From this, we see that if assumption 
(ii) holds then AMS measurements should give consistently older radiocarbon dates than those 
obtained using radiometric methods for > 3 cal kBP. 

While the predicted A 14 C curve which arises from case (i) cannot easily be verified experi- 
mentally, our prediction from case (ii), that radiocarbon dates obtained using the two different 
methods should be significantly different for older samples, is very easy to check using results 
of the Fourth International Radiocarbon Intercomparison (FIRI) [21]. These results show no 
significant difference in the predicted ages of a range of standard samples, including a sample 
of humic acid (sample E) dated at around 11.8 kBP, which is just beyond the range that we are 
considering. The AMS mean age for this sample was 54 ±53 14 C yr higher than the correspond- 
ing age obtained using Gas Proportional Counting, a radiometric method, which is nowhere 
near the approximately 3,000 14 C yr difference that our results predict for this sample. 

These results would appear to invalidate the assumptions of case (ii) and support the con- 
ventional (but not verified) assumptions of case (i). However, before dismissing case (ii), we 
make two proposals that might explain why this predicted difference in dates is not observed. 

Firstly, case (ii) is the extreme case. In the more likely scenario in which the discrepancy 
in the calibration curve is due to a combination of atmospheric variation and non-exponential 
decay the predicted difference between Cams and Cp would be smaller. 

Secondly, we note that radiometric methods were the first to be used for carbon dating. 
In the late 1970's, AMS techniques were developed to measure the 14 C/ 12 C ratio directly [25J. 
However, 'absolute AMS isotopic measurements are extremely difficult' [26J and so these meth- 
ods were calibrated using international standards and modern samples of known age. For older 
samples, and given a strong belief in exponential decay, the results were also compared with 
the predicted radiocarbon ages obtained by radiometric methods. Indeed, the first aim of the 
FIRI exercise was 'demonstration of the comparability of routine analyses carried out by both 
AMS and radiometric laboratories' [2"3] . 

In practice, AMS measurements are the product of complex experimental procedures [25J. 
Results are affected by a number of factors, such as contamination during chemical processing 
[7J [271 [28], sample size [27J [29] and pretreatment of samples by acid hydrolysis [291 121] • A 
small change in the AMS measurement as a result of any of these factors would be sufficient to 
generate a significant difference in the predicted radiocarbon age for older samples. 

We therefore conclude that as long as AMS measurements are calibrated to give agreement 
with results from radiometric methods, our predicted difference in radiocarbon dates between 
the two methods will never be observed. 



T r = CAMs(Td) 
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8 Conclusions 



The question posed in the title of this paper as to whether radioactive decay is really exponen- 
tial (in all circumstances) is of fundamental importance, both theoretically and practically. We 
have presented experimental evidence, from the results of radiocarbon dating and modelling, 
which suggests that non-exponential decay may be a significant factor in the discrepancy in 
radiocarbon dates. We have also noted that quantum mechanics does not predict exact expo- 
nential decay, and have presented mathematical arguments that non-exponential decay should 
be expected for slowly decaying isotopes, but not for rapidly decaying isotopes. Thus, we be- 
lieve that there is strong evidence, both experimental and theoretical, that radioactive decay 
of slowly decaying isotopes is not exactly exponential. 

We have also shown that one significant consequence of non-exponential decay is that the 
two methods used for radiocarbon dating, namely radiometric and AMS, should give different 
results for older samples. This effect is not seen in practice, but we suggest that this could 
be due to AMS being calibrated to give agreement with results obtained using radiometric 
methods. 

Our results in the previous section provide a way to experimentally test our prediction that 
the discrepancy in the calibration curve is due to a combination of atmospheric variation and 
non-exponential decay. To do this, a careful calibration of AMS measurements is required to 
ensure that correct absolute measurements are being made, particularly in the range 3-12.55 
cal kBP. Dating of tree samples in this range should then be performed. If these radiocarbon 
dates are significantly higher than those in the IntCal09 dataset, then it would confirm our 
prediction. Conversely, if good agreement is obtained with the IntCal09 data, then it would 
disprove our prediction. 

The question we are addressing of whether radioactive decay is exponential over long time 
periods is certainly of much theoretical interest. However, the stakes are high in more practical 
realms, since many dating methods over long time periods rely on the assumption of exponential 
decay of a slowly decaying isotope [I]. If the decay of 14 C is indeed non-exponential, then 
no other slowly decaying isotope can be assumed to decay exponentially either, which would 
remove a foundation stone of modern dating methods. If confirmed, this will require a radical 
reappraisal both of our approach to radioisotope dating methods and of the many currently 
accepted dates which have been obtained using these methods. 
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